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abstract 

We suggest that a neutron star with a strong magnetic field, spiraling into the envelope of a companion 
star, can generate a "companion induced SN-like off-center explosion" . The strongly magnetized neutron 
star ("magnetar") is born in a supernova explosion before entering into an expanding envelope of a 
supergiant companion. If the neutron star collapses into a black hole via the hypercritical accretion 
during the spiral-in phase, a rapidly rotating black hole with a strong magnetic field at the horizon 
results. The Blandford-Znajck power is sufficient to power a supernova-like event with the center of 
explosion displaced from the companion core. The companion core, after explosion, evolves into a C/O- 
white dwarf or a neutron star with a second explosion. The detection of highly eccentric black-hole, 
C/O-white dwarf binaries or the double explosion structures in the supernova remnants could be an 
evidence of the proposed scenario. 

Subject headings: binaries: close - stars: neutron - stars: evolution - stars: statistics 



1. INTRODUCTION 

It has recently been reported that neutron stars with 
large luminosities could be powered by ultra-strong mag- 
netic fields ~ 10 14 G ("magnetars") rather than rapid ro- 
tation as in rotation-powered pulsars (e.g. Kouveliotou 
et al. 1998, 1999). These "magnetars" are believed to 
be produced in supernova explosions or accretion induced 
collapses (Duncan & Thompson 1992). 

On the other hand, the Blandford-Znajek process has 
been suggested as a powerful source for a variety of very en- 
ergetic events ranging from active galactic nuclei (Bland- 
ford & Znajek 1977) to gamma-ray bursts (e.g. Lee et al. 
2000ab and references therein). The extremely high lu- 
minosities or power outputs generally require ultra-strong 
magnetic fields around rapidly rotating black holes. The 
origin of such strong magnetic fields around the rapidly 
spinning black holes has not been clearly understood. It 
is interesting to point out that if magnetars directly col- 
lapse to black holes via rapid accretion, the resulting black 
holes could spin rapidly with magnetic fields left at the 
horizon, which are large enough to make the Blandford- 
Znajek power interesting for very energetic events. 

Based on this possibility, we propose that the hyper- 
critical accretion onto a neutron star with a very strong 
magnetic field can produce a rapidly spinning black hole 
with an ultra-strong field near the horizon. The hyper- 
critical accretion has been invoked for the common enve- 
lope phase evolution of the massive stellar binaries (Bethe 
& Brown 1998). Unless the strong magnetic fields decay 
within time scales shorter than the duration of the accre- 
tion phase, ~ 1 year, it is highly likely that the black holes 
with very strong fields could be produced. The Blandford- 
Znajek power is sufficient to power a supernova-like event 
with the center of explosion displaced from the compan- 
ion core. The companion core, after explosion, evolves 
into a C/O-white dwarf or a neutron star with a second 
explosion. The detection of highly eccentric black-hole, 



C/O-white dwarf binaries or the double explosion struc- 
ture in the supernova remnants, with two different centers 
of explosion, could be an evidence of the proposed sce- 
nario. Especially, in the later case, the neutron star can 
have double kicks, resulting in very high kick velocities, 
> 1000 km s _1 . If the spiral-in of the neutron star leads 
to the merger with the companion core, the explosion at 
the core would appear as an explosion from a single star. 
The neutron star-neutron star binary system is expected 
when the companion collapses before significant accretion 
of the neutron star occurs. 

2. MAGNETAR FORMATION AND ULTRA-STRONG 
MAGNETIC FIELDS 

The NSs right after core collapse are hot and convective 
throughout most of the stellar interior due to a large neu- 
trino flux. Duncan & Thompson (1992) and Thompson & 
Duncan (1993) describe the most likely outcome of such 
hot neutron stars. Right after the collapse, there exists 
an intense neutrino-driven convective phase with the over- 
turn time scale of convective cell, r conv ~ 1 ms at the base 
of the convection zone and r conv ~ 0.1 ms at the neutri- 
nospherc. In this case, if the neutron star spin frequency 
is 
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the amplification of of a magnetic field by helical motions 
is not suppressed by turbulent diffusion, and the aft-type 
dynamo could work efficiently. They found that the sat- 
uration field strength of £? sat ~ 10 16 — 10 17 G could be 
reached within ~ 10 — 100 ms. Therefore, post-collapse 
NSs with P 1ms, within the first few seconds, can build 
up very strong large scale magnetic fields > 10 15 G through 
exponential growth. In this case, the resulting ultra-strong 
field is independent of the fossil field strength. For strong 
seed fields, the small-scale fluid motion, which is required 
for ohmic dissipation of magnetic energy, is suppressed by 
the magnetic tension. The convective motion ceases when 
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the neutrino-driven convection becomes unavailable. Such 
a strong field pulsar with slow spin (" magnetar" , Duncan 
& Thompson 1992) has recently been discovered, which 
strongly supports the basic idea of creation of a ultrastrong 
field pulsar (Kouvelioutou et al. 1998, 1999). 

For NSs spinning with Q <^ ^dynamo or P > 1ms, the dif- 
ferential rotation (i.e. dQ/dR ^ 0) could amplify the seed 
poloidal fields by wrapping them up into strong toroidal 
fields (Kluzniak & Ruderman 1998). For an initial poloidal 
seed field with strength B p ~ 10 9 G which would result 
from flux freezing of a field of 3 x 10 4 G in the pre- 
collapse core, is amplified to the strong toroidal field of 
P - 10 17 G after - (l/27r)(P /P p ) - 10 7 rotations or 
time t ~ B^/BpVl ~ 10 4 P s with P in ms. The ampli- 
fied field becomes buoyant after reaching the critical field 
strength of B\, ~ 10 17 G. Since the total energy included 
in the buoyant fluid cells is - ^-B\{^-)t] ~3x 10 51 ?7 
erg, where rj ~ 0.2 is the fractional volume of the buoy- 
ant elements, the repeated build-up and buoyant loss will 
occur for a number of times determined as 
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where / s=s |MP 2 with neutron star mass M = 1.5 Mq 
and neutron star radius R = 10 km. This process will con- 
tinue until the existing poloidal seed field is exhausted as 
the poloidal field is not re-generated in this non-dynamo 
amplification process. The buoyantly rising field could dis- 
sipate on a time scale 
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which controls the overall time scale. Therefore, we con- 
clude that essentially all NSs with P ~ a few ms could 
achieve strong fields through either dynamo or differential 
rotation. 

Once strong fields are created, rapid spin-down of NSs 
is inevitable. The electromagnetic dipole radiation lumi- 
nosity is 
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which gives the characteristic electromagnetic spin-down 
time scale 
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where B 15 = P/10 15 G and Q t = ft/10 s- Therefore, 
all MSPs with ms periods are likely to spin-down to near 

1 /2 

the pulsar deathlinc, P ~ 705,5 s, on a time scale of 
~ 10 5 B^ yr (e.g. Ritchings 1976). That is, these initially 
fast spinning MSPs would not be detected as MSPs unless 
they are recycled through the low mass X-ray binary stage, 
which is not promising (Yi & Grindlay 1998). Therefore, 
before the common envelope phase begins, the neutron 
stars are likely to be spinning with spin periods very simi- 
lar to those of the magnetars and anomalous X-ray pulsars 



(e.g. Mereghetti 1999, Kouvelioutou 1998, 1999). That is, 
the initial conditions for the neutron stars which spiral-in 
during the common envelope phase are likely to be a slow 
spin and a high magnetic field. 

3. HYPERCRITICAL ACCRETION 

Bethe & Brown (1998) suggested that the hypercriti- 
cal accretion can occur in the massive binary systems, in 
which both stars can go into supernova leaving the com- 
pact cores as remnants. In their scenario, the more massive 
star 1 first explodes leaving a neutron in the core. If the 
binary is within the radius of the companion in the giant 
stage, the first born neutron star will spiral into the enve- 
lope of the expanding companion in the giant stage. We 
assume the envelope of the giant to be convective. 

In the rest frame of the compact object, Bondi-Hoyle- 
Lyttleton accretion of the hydrogen envelope of density 
Poo and velocity V is 
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where M co ,i = M co /l Mq is the mass of the compact ob- 
ject, Vg is the velocity in unit of 1000 km s^ 1 , and is 
given in g cm -3 . The density poo required for the hyper- 
critical accretion, M = l0 s M Edd w 10 26 g s _1 , is 



V 3 

Poo - 0.44 x Kr 3 -rrf- g cm 



-3 



Ml,, 



(7) 



Using the Kepler law for circular orbits V 2 = GM tot /a 7 
where M tot is the combined mass of the compact object 
and the stellar material interior to the orbit of the compact 
object and a is the binary separation, one finds 
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where M tot ,io = M tot /10 10 M Q and a i2 = a/10 12 cm. Since 
this critical density is of the same order as the average 
density of the hydrogen envelope of the giant star, the 
hypercritical accretion M ^ lO 8 -/^^ can occur in this 
spiral-in process. This rapid accretion can add ~ 1 Mq 
to the first-born neutron star in less than about one year, 
which is enough to turn the first born neutron star into a 
black hole. In this process, the orbital radius shrinks by a 
factor ~ 50. 

If the companion is massive enough to have a super- 
nova explosion, the binary system will end up as a black 
hole— neutron star (bh,ns) binary instead of a (ns,ns) bi- 
nary. If the companion is less massive, the binary will 
end up as black hole— CO white dwarf (bh,cowd) bina- 
ries. Brown, Lee, Portegies Zwart, & Bethe (1999) show 
that there are big discrepancies between the predicted 
and the observed (ns,ns) and (ns,cowd) c (i.e. with circu- 
lar orbits) binary systems. They solved this discrepancy 
by arguing that in both systems the first born neutron 
star turned into black holes by accretion leaving behind 



1 The progenitor masses should be different by more than 5%. If the two masses are close enough, the binary will explode nearly at the 
same time leaving the neutron star— neutron star (ns,ns) binary. In this case, there is not enough time to change neutron star into black hole 
by accretion. 



3 



(bh,ns) and (bh,cowd) binaries. The (bh,ns) binaries are 
one of the most important sources for gravitational detec- 
tors. The chirp mass of double neutron star binary with 
masses 1.4 M© is M chirp = ^ S {M 1 + M 2 ) 0A = 1.2 M , 
while that of the (bh,ns) binary with Mbh = 2.4 Mq is 
M c hirp = 1-6 Mq. The (bh,ns) binaries, if they exist, can 
be more easily detected by factor (1.6/1.2) 3 = 2.4 than the 
(ns,ns) binaries. Since the LIGO is based on one observa- 
tion of (ns,ns) binary, their unseen (bh,ns) will enhance 
the detectability of LIGO by factor of ^30 considering the 
effect of the chirp masses. This will be tested in near fu- 
ture. 

One of the most strong criticisms on the hypercritical 
accretion was the angular momentum and possibility of a 
strong jet along the rotation axis. Even thouth these re- 
main as unsolved problems, Brown, Lee, & Bethe (2000) 
showed that Narayan & Yi (1994) ADAF solution can be 
applied in this case if the accretion flows is radiation- 
dominated. With the hypercritical accretion, the in- 
ner boundary of accretion flow can be extended to the 
marginally bound orbit (r m & = 2Rs c h for Schwarzschild 
black holes). Since r m & « 9 km, the inner boundary of 
accretion disc can be extended to the neutron star surface 
fNS ^ 10 km. By using a reasonable viscosity parameter 
a = 0.05, they found that the neutrino cooling in the disc 
is almost negligible < 10~ 7 of viscosity generated energy. 
Brown ct al. (2000) argued that in order to avoid strong 
polar jets or pile-up up of material above neutron star sur- 
face, the neutrinos are generated at the neutron surface to 
carry out the pressure and energy from the system. 

4. HYPERCRITICAL ACCRETION ONTO A MAGNETAR 

Despite a very strong magnetic field, the dynamic ef- 
fects of the magnetic field would be negligible if the mass 
accretion rate satisfies 

(M \ 1//2 

~ 4 x 10 7 B 2 14 [-^rj M Edd (9) 

Such a possibility is realized when the accretion is hyper- 
critical. 

During hypercritical accretion, even when the electro- 
magnetic dipole radiation occurs for a fast spinning pul- 
sar, the spin-up by hypercritical accretion can counter the 
spin-down. The spin-down torque due to the dipole radi- 
ation emission is 

N em = 1-B 2 R 6 sin 2 d fi 3 (10) 
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and the spin-up torque due to accretion is 

N acc = Ma {GMR) 1 ' 2 (11) 

where a ~ 0.5 is the sub-Keplerian rotation factor in the 
hypercritical accretion (Brown, Lee, & Bethe 2000). The 
equilibrium spin is reached at the neutron star spin fre- 
quency 
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By assuming (sin 2 9) ~ 0.5, one gets VLns ~ 1-8 x 10 3 s _1 
with other quantities in unities. The collapse of the neu- 
tron star to a black hole results in a rapidly spinning black 
hole with a spin frequency 

which corresponds to black hole spin of a ~ 0.1. The 
magnetic field strength 

Bh = (jl^j B m » 4B NS w 4 x 10 14 B NS , 14 G (14) 

where we have assumed the angular momentum conserva- 
tion and flux-freezing. 

If the magnetar enters the hypercritical accretion regime 
as a slowly spinning neutron star, the spin-up during the 
hypercritical accretion will take place on a time scale 

tspin-up 1.6 x 10' s - 0.5 yr (15) 
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which is comparable to the spiral-in time scale (~ 1 yr) 
during the common envelope phase. 

5. ENERGETICS OF BLANDFORD-ZNAJEK PROCESS 

Based on the steady pulsar periods and spin-down rates, 
the magnetic field of pulsars are generally estimated by 
assuming the dipole radiation induced by magnetic fields. 
The typical fresh neutron stars in the pulsar island have 
~ 10 12 G, and the recycled ones have <^ 10 10 G. The dis- 
coveries of the magnetar and AXP suggest the possibilities 
of the strong magnetic field J> 10 14 G. Even though AXP 
and magnetars raised questions on the origin of the pul- 
sar mechanism, we assume the dipole radiation as origin 
of the pulsar in these systems. The probability of pulsar 
with strong magnetic field is assumed to be ~ 10% (Kou- 
veliotou et al. 1998, 1999). 

After the discovery of GRB afterglows, the theoretical 
understandings of the formation of GRBs have been im- 
proved quite a lot. However, the mechanism of the central 
engine is still controversial. The electromagnetic extrac- 
tion of the BH rotating energy by the Blandford-Znajek 
process (Bandford & Znajek 1977) is one of the candidate 
mechanisms, in which the strong magnetic field is essen- 
tial to power the GRB. Typically Bh ~ 10 15 G can gener- 
ate enough energy to power the GRB from rotating black 
holes. However, the existence of the strong magnetic field 
in the black holes is uncertain because the strong magnetic 
fields diffuse in very short time scales. 

We suggest that the Blandford-Znajek process can oper- 
ate if the collapsing neutron stars, with B 5> 10 G, carry 
the magnetic fields with them, amplifying the fields upto 
~ 10 15 G. In this scenario, the existence of neutron stars 
prior to the black holes is essential. 

The hypercritical accretion prior to the collapse to a 
black hole spins up the neutron star and the black hole 
is expected to be spinning at rotational angular velocity 
close to a ~ 0.1. Once black hole is formed, the spin-down 
torque due to the dipole emission is not working, and the 
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spin-up torque due to accretion can increase the black hole 
spin further. So, our estimate given above may be a low 
limit. 

The Blandford-Znajek power (Lee, Wijers, & Brown 
2000a) from the rotating black hole with the angular mo- 
mentum J = aM and the magnetic field Bh is given by 

Pbz ~ 10 50 a 2 B 2 H15 (M/M Q ) 2 erg s~\ (16) 

Using the values in the previous section, we have Pbz ~ 
3.6 x 10 47 erg s _1 . The total available rotational energy 
is 0.06% of the rest mass ~2x 10 51 erg (Lee, Wijers, & 
Brown 2000a). 

6. THE ROLE OF BLACK HOLE FORMATION IN THE 
SPIRAL-IN PHASE 

When the first born star with strong magnetic field ac- 
cretes material from the companion envelope, the black 
hole should be formed inside the envelope. Since the black 
hole is formed in the process of the hypercritical accre- 
tion, the temperature near the surface of the neutron star 
should remain ~ 1 MeV until the collapse into a black 
hole. This temperature is required to generate neutrinos. 
In this case, the spin frequency for the r-mode instability is 
so high that we can neglect the effect (Lindblom & Owen 
1999). If the Blandford-Znajek process works at the time 
of black hole formation, there is enough energy to generate 
GRBs. However, because of the optically thick envelope, 
gamma rays cannot come out through the hydrogen (or 
helium) envelope. Instead of GRBs, the energy will pile 
up in the envelope to generate supernova. The main dif- 
ference from the normal SN scenarios is that the SN is 
not induced by the core collapse of the progenitors. The 
natural consequence is the strong asymmetry in the SN 
events. Since the companion mass is not the unique pa- 
rameter, even the less massive stars (normal white-dwarf 
progenitors) can have SN events with high kick velocities. 

The available energies powered by rotating black holes 
(> 1.5 Mq) with Blandford-Znajek process can be ~ 10 51 
erg. Since the binding energies (O.QGM 2 / R) of the hydro- 
gen envelope of giant stars (10 M ) are ~ O(10 48 ) erg, 
newly formed black holes have enough energy to blow the 
envelope off. If only 1% of this SN energy is available for 
the kinetic energy of the He core of the giant star, the 
available kick velocities are ~ O(1000 km s _1 ). 

Fast moving white dwarfs from SN remnants, if exists, 
will be the consequences of this "companion induced SN 
events" . If the black holes are formed just before the nor- 
mal SN events of the companions, the double SN are also 
possible leaving (bh,ns) as remnants. Especially, in the 
later case, the neutron star can have double kicks, result- 
ing in very high kick velocities, > 1000 km s _1 . 

Bcthe & Brown (1998) argued that the formation rate 
of (bh,ns) binaries is 10~ 4 yr _1 per galaxy, where the low 
mass black holes (~ 2.5 M©) are formed in the common 
envelope evolution. By assuming the flat q distribution, 
Brown, Lee, Portegies Zwart, & Bethe (1999) also argued 
that the formation rate of circular (bh, cowd) c binaries is 
~ 1.7 x 10~ 4 yr^ 1 per galaxy. In their scenario the first 
born neutron star went into the common envelope evolu- 
tion, where the neutron star is converted into black hole 
by hypercritical accretion. According to our scenario, in 
the formation of (bh,ns) and (bh,cowd) c binaries by hy- 
percritical accretion, about 10% of the neutron stars have 



very strong magnetic field and will have the companion 
(newly formed black holes) induced supernova-like explo- 
sions. As a result ~ 10% of progenitors of (bh,cowd) c 
binaries end up as eccentric (bh, cowd) e binaries instead 
of circular ones. Also ~ 10% of progenitors of (bh, ns) 
binaries have double explosions with two different centers 
of explosion. By assuming that half of the progenitors 
of (bh, cowd) e or (bh, ns) binaries survive the explosion in- 
duced by rotating black holes, we have the formation rates 
of eccentric (bh,cowd) e binaries or double explosions 

R ~ 10~ 5 yr _1 per galaxy. (17) 

Since the cooling time of CO white dwarfs, i.e. the time re- 
quired to reach luminosities log(L/L©) = —4.5 with Lq = 
solar luminosity, is 1 — 10 Gyr (Salaris et al. 1997), the 
number of eccentric (bh,cowd) e binaries with luminosity 
log(L/L©) > -4.5 is 

N ~ 10 4 - 10 5 in Galaxy. (18) 

The number of non-pulsating eccentric (ns, cowd) e bina- 
ries with luminosity log(L / Lq) > —4.5 is ~ 10 times more 
popular than (bh,cowd) e systems. In order to identify 
(bh,cowd) e systems, one need well-established maximum 
neutron star masses. If there are mechanisms to identify 
the black holes from dead pulsars in addition to the mass 
(Balbcrg ct al. 1999), one can detect (bh,cowd) e systems. 

7. DISCUSSION 

We have looked into a possibility that a supernova-like 
explosion could occur during the common envelope evo- 
lution of a strongly magnetized neutron star. The "mag- 
netar" accretes in the form of the hypercritical accretion 
during the spiral-in phase and collapses to a rapidly spin- 
ning black hole with a strong magnetic field near the hori- 
zon. Then the explosion results from the sudden release 
of energy through the Blandford-Znajek mechanism. In 
this scenario, the strong magnetic field naturally occurs 
as it is derived from the strong field of the pre-collapse 
" magnetar" . 

There are a number of points to be clarified in the pro- 
posed scenario some of which could be observationally 
testable. First, if the magnetic field of the magnetar de- 
cays on a time scale ~ 10 3 — 10 4 yrs as is often argued 
for the evolutionary behavior of magnetars and anoma- 
lous X-ray pulsars (Colpi, Geppert, & Page 2000), then 
the common envelope phase has to follow the magnetar 
formation. Second, during the accretion phase, most of 
the emitted energy (including the dipole radiation) is ul- 
timately radiated from the neutron star in the form of the 
neutrinos. Third, if the Blandford-Znajek phase occurs 
while the neutron star collapses into black hole in the en- 
velope, the resulting off-center explosion will give a high 
eccentricity, which should be testable in terms of the highly 
eccentric binary systems. Fourth, there will be little metal 
spread from the core since there will be no core collapse. 
The dominant effect would be blowing-off the hydrogen 
envelope preceded by the expansion of the envelope. 

If the merger of the companion core with the neutron 
star occurs before the black hole formation, then the ex- 
plosion will not cause an off-center explosion. Although 
this possibility exists, the estimated accretion time scale 
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on which the neutron star can accrete enough mass to col- 
lapse to a black hole is comparable to the spiral-in time 
scale. This implies that the off-center explosion is possible. 
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